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Abstract 

Competing pinning effects on a D- dimensional interface by weak impurity 
disorder and a periodic potential of the underlying crystal lattice are 
analyzed for 2 < < 4. We use both the Gaussian variational method 
(GVM) and the functional renormalization group e = A — D expansion 
(FRO) which yield different phase diagrams: Whereas the FRG always 
predicts a rough phase with irrelevant lattice pinning, the GVM in 
combination with a three parameter RG for the random potential [T. 
Nattermann, H. Leschhorn, Europhys. Lett. 16 (1991) 603] leads to 
a roughening transition of first order. For random bond disorder we 
compute self-consistently the effective lattice potential. 
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Einleitung 



Oberflachen und Grenzflachen spielen in der gesammten (naturwissenschaftlichen) Welt 
eine groBe RoUe. Beispiele hierzu gibt es aus der Biologie, der Chemie und, was uns 
besonders interessiert, aus der Physik. Der wichtigste Vertreter aus der Biologie ist die 
Zellmembran, die es iiberliaupt ermoglicht, das es Lcbcn gibt. Aus der Chemie ist sicherlich 
der Katalysator ein wichtiges Beispiel, den wir allc aus dcm Auto kennen. Er bcschleunigt 
chemische Reaktionen, indem er, durch Anlagerung der Schadstoffe an der Oberflache, es 
erleichtert, deren chemische Bindungen aufzubrechen. So konnen neue, weniger schadliche 
Stoffe enstehen. Sein Wirkungsgrad ist proportional zur effektiven OberflachengroBe. 

Aus der physikahschen Sicht ist eine Oberflache ein Objekt, welches zwei Telle eines 
Systems voneinander separicrt. Die Dimension der Oberflache ist immer um eins kleiner 
als die des zugrundeliegenden Raumes. Die Mathematiker sprechen in einem solchen Fall 
von einer Hyperflache. Beispiele gibt es ebenso unterschiedliche, wie zahlreiche: So gibt es 
"gewohnliche" Oberflachen wie die Begrenzungen eines Festkorpers. Die Grenze zwischen 
zwei koexistierenden Phasen ist ebenfalls eine Oberflache. 

Die Physik an der Oberflache unterscheidet sich vielfach von der in dem dazugehorigen 
Volumcn. Bci Festkorpern kommt es vor, dafi das Gitter an der Oberflache eincn andcrc 
Struktur hat. Im einfachsten Fall hat die oberste Atomlage einen anderen Gittcrabstand 
senkrecht zur ihrer Ausbreitungsrichtung, als die Lagen tief im inneren des Kristalls. Im 
Extremfall besitzt die Oberflache eine veranderte Gitterstruktur, indem sie eine Uberstruktur 
ausbildet. Man spricht von einer Rekonstruktion der Oberflache, bei der die Gitterkonstante 
parallel zur Oberflache ein Vielfaches der Gitterkonstante im Volumcn ist. Zudcm kann 
dieses in verschienen Richtungen ein unterschiedliches Vielfaches sein. Die Bandstruktur 
der Oberflache unterscheidet sich ebenfalls von der im Volumen. Somit andern sich auch 
die elektronischen Eigenschaften. 

Die modernen experimentellen Methoden ermoglichen es heute eine atomare Oberflache 
direkt zu beobachten. Speziell erwahnt seien hier das Raster-Tunnelelektronen-Mikroskop 
und das Raster-Kraft-Mikroskop. Erstercs crlaubt cs die elektronische Struktur der Oberflache 
mithilfc eines Tunnelstroms, der durch cine winzigc, im idealen Fall einatomige Spitze flieBt, 
abzutastcn. Das Raster-Kraft-Mikroskop tastct die Oberflache direkt mit einer Spitze ab 
und gibt somit Aufschlufi iibcr den atomarcn Aufl^au. 

In dieser Arbeit geht es darum einen der vielen Aspekte, welche eine Oberflache besitzt, 
etwas genauer zu beleuchten. Wir woUen deren Rauhigkeit betrachten. Es handelt sich 
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dabei um eine Gleichgewichtseigenschaft. Diese Eigenschaft hat auch ihre Auswirkungen 
auf Nichtgleichgewichtsphanomene. Als Beispiele seien hier getriebene Grenzflachen, das 
Aufdampfen, auch bekannt unter dem enghschen Namen Molecular Beam Epitaxy (MBE), 
oder allgemeiner das Kristallwachstum genannt. 

Das hier betrachtete Modell ist nicht nur als Oberflachenmodell interessant. Das reine 
Modell gehort zu einer ganzen Gruppe von 2-dimensionalen Modellen, die sich aufeinander 
abbilden lassen, wie das XY — Modell und das Coulombgas. Weiterhin gibt es mit diesem 
Modell verwandte Probleme. Die Oberflache beschreiben wir durch ein skalares Feld, welchen 
durch einen in der Ebene liegenden Vektor beschrieben wird. Andern wir die Dimension 
dieses Parametervektors z.B. auf den Wert eins und verwenden wir ein zweikomponentiges 
Felds, so erhalten wir die Beschreibung eines gerichteten Polymeres. Mit einem ganzen 
Ensemble solcher Linien kann man das FluBliniengitter eines Typ 11 Supraleiters in der 
Abrikosov-Phase beschreiben. Als eine andere Variation kann man ein Ensemble von Grenzflachen 
betrachten, womit wir bei der Theorie der kommensurablen Systeme angekommen sind. 

Die Arbeit ghedert sich folgendermaBen: Zunachst wird ein Modell motiviert, mit dem 
man Oberflachen beschreiben kann. Es wird auf einige bekannte Resultate von verwandten 
Modellen eingegangen. Im darauffolgenden Kapitel wird das Modell mithilfe der Renormierungsgruppe 
untersucht. Im Rahmen einer e-Entwichlungen werden FluBgleichungen gewonnen, aus 
denen wir erste Resultate gewinnen. Danach folgt eine Behandlung mittels eines Variationsverfahrens. 
Hierbei wird eine Abwandlung der Renormierungsgruppe aus dem vorangehenden Kapitel 
verwendet, die Drei- Parameter- Approximation. SchlieBlich folgt ein Kapitel, in dem eine 
Methode vorgestellt wird, mit der man auch in der urspriinglichen Renormierungsgruppe 
zu einem Phaseniibergang gelangen kann. 

Fiir die Beschreibung einer Oberflache hat man zwei Moglichkeiten. Die eine besteht darin, 
alles auf einer atomaren Langenskala zu betrachten. Hierbei wird die Position eines jeden 
Atoms angegeben. Ein bekanntes Beispiel fiir diese Art der Beschreibung sind ist das Solid 
on Solid (SOS) Modell. Die andere Moglichkeit ist, iiber kleine Bereiche zu mitteln und 
somit zu einer Kontinuumsbeschreibung iiberzugehen. Der zweite Zugang wird in dieser 
Arbeit angewendet. 
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Kapitel 1 
Das Modell 



1.1 Ein Oberflachenmodell am Beispiel des Ising-Modells 

Als Ausgangspunkt fiir unsere Betrachtungen wahlen wir das Ising-Modell auf einem 
Quadratgitter. Wir stellen uns vor, daB auf jedem Gitterpunkt ein Spin mit zwei Einstellmoglichkeiten 
sitzt, der Einfachheit halber s = ±1. Diese Spins konnen nun jeweils mit ihren nachsten 
Nachbarn wechselwirken. Die Starke der Wechselwirkung wird durch eine Kopplungskonstante 
J beschrieben. Die Hamiltonfunktion fiir ein reines System, das heiBt ein System, wo die 
Kopplungskonstanten zwischen den nachsten Nachbarn alle gleich grofi sind, hat damit die 
folgenden Gestalt 

"Rising ^ ^ Js^Sj. 

<i,j> 

< i,j > deutet an, daB nur iiber nachste Nachbarn summiert wird, wobei jedes Paar nur 
einmal gezahlt wird. (Anderenfalls miiBte man den Faktor 1/2 erganzen.) 

Fiir dieses System gibt es nun zwei Falle mit je zwei entarteten Grundzustandskonfigurationen, 
d.h. Konfigurationen mit minimaler Energie E. Die T = Konfigurationen sind diese 
Grundzustandskonfigurationen. 1st die Kopplungskonstante J > 0, so richten sich die 
Spins parallel zueinander aus. Dieses ist der ferromagnetische Fall. Der Grundzustand ist 
zweifach entartet, da das System noch die Moglichkeit hat, alle Spins nach oben oder nach 
unten zeigen zu lassen. Bei J < 0, dem antiferromagnetischen Fall, sind benachbarte Spins 
entgegengesetzt ausgerichtet. Das Gitter besteht aus zwei Untergittern, in denen jeweils die 
Spins parallel ausgerichtet sind. Auch hier gibt es zwei mogliche Konfigurationen, entweder 
zeigen in dem einen Untergitter die Spins nach oben und in dem anderen nach unten oder 
eben umgekehrt. Wir betrachten im folgenden den Ferromagneten. 

Vom Ising-Modell ist bekannt, daB es in zwei und mehr Dimensionen einen Phaseniibergang 



besitzt |T^. Unterhalb der kritischen Temperatur ist die Magnetisierung, das ist der 



thermisch gemittelte Erwartungswert des Spins, ungleich null, oberhalb gleich null. Der 
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Phaseniibergang ist von zweiter Ordnung, das heiBt die Magnetisierung verschwindet kontinuierlich. 
Starten wir zum Beispiel bei T = und erhohen die Temperatur, so findet ein in der 
statistischen Physik immer wieder auftretender ProzeB statt. Das Gleichgewicht ist durch 
das absolute Minimum der Freien Energie T = E—TS bestimmt {S ist die Entropie). Einige 
Spins drehen sich um, wodurch sich die Energie des Systems erhoht. Auf der anderen Seite 
gewinnt das System aber an Entropie, da diese Spins an beliebigen Platzen im System liegen 
konnen. Mit wachsender Temperatur bilden sich aus vereinzelten 'verkehrt' stehenden Spins 
kleine Bereiche, sogenannte Domanen. Der typische Durchmesser dieser Domanen, welchen 
man als die Korrelationsldnge bezeichnet, wachst, bis er bei die SystemgroBe erreicht. 
Genauer muB man sagen, da es einen Phaseniibergang nur im thermodynamischen Limes 
gibt, daB die Korrelationslange bei Tc divergiert. Unter dem thermodynamischen Limes 
versteht man, daB man die SystemgroBe gegen Unendhch gehen laBt, wobei man die Dichten 
des Systems konstant halt. 

Was uns nun interessiert sind die Begrenzungen dieser Domanen, die sogenannten Doma- 
nenwande. Fiir sie wollen wir eine geeignete Hamiltonfunktion aufstellen und damit deren 
Eigenschaften untersuchen. Zunachst machen wir uns ein paar Gedanken iiber die Energie 
einer solchen Wand. Die Wand befindet sich genau da, wo die Spins entgegengesetzt stehen. 
Im Vergleich zu parallel ausgerichteten Spins ist die Energie an solchen Stellen um 2 ■ J 
groBer. Die Anzahl solcher gebrochener Bindungen ist proportional zur Flache dieser Wand. 
Somit setzen wir die Energie der Wand als eine Flachenenergiedichte multipliziert mit der 
Flache der Wand an. 

Was wir noch benotigen ist eine geeignete Parametrisierung der Domanenwande. In dem 
oben beschriebenen System ist dieses mathematisch sehr kompliziert. Wir wollen uns daher 
auf eine einzige Domanenwand beschranken. Eine solche Wand kann man im Ising-Modell 
erzwingen, in dem man in einer Raumrichtung antiperiodische Randbedingungen fordert. 
Diese ausgezeichnete Richtung wollen wir im folgenden als die z-Richtung bezeichnen. 
Durch diese Randbedingung muB das System im Verlauf der z-Richtung die Spinrichtung 
wechseln. Im Grundzustand wird eine vollig ebene Wand vorliegen. Eine solche Wand 
laBt sich nun einfach durch orthogonale Koordinaten parametrisieren (Abbildung |1.1| ). Die 
D := d—1 Koordinaten, entlang der Achsen senktrecht zur Richtung bezeichnen wir mit 
dem Vektor x, die entlang der z-Achse mit z. Die zu untersuchende Flache laBt sich nun als 
Feld z{x) beschreiben. Die Punkte der Flache haben demnach die Koordinaten {x,z{x)). 
Hierbei bezeichnen wir mit d die Dimension des ganzen Raumes und mit D die Dimension 
der Flache. Eine Flache, welche so parametrisiert ist, bezeichnet man als eine gerichtete 
Flache. 

Den Preis, welchen man fiir diese einfache Parametrisierung zu zahlen hat, ist, daB man mit 
ihr keine Uberhange oder Einschliisse beschreiben kann. Die Frage ist nun, ob dieses unsere 
Betrachtungen wesentlich beeinfiuBt. Wir gehen zu einer vergroberten Beschreibung (coarse 
graining) des Systems iiber, das heiBt, wir fassen einen ausgedehnten Bereich der Oberfiache 
zusammen und summieren iiber die Fluktuationen in diesem Bereich. Wenn wir diesen 
Bereich groB genug wahlen, so werden die Uberhange und Einschliisse herausgemittelt. 
Wir erreichen hierdurch eine Kontinuumsbeschreibung der Oberfiache. Der EinfiuB der 
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1.1 Ein Oberflachenmodell am Beispiel des Ising-Modells 




Abbildung 1.1: Die Koordinaten einer gerichteten Fldche: der Parametervektor x spannt 
eine D-dimensionale Hyperebene auf; das Feld z{x) beschreibt das Hdhenprofil der 
Oberfldche iiber dieser Hyperebene 

Einschliisse und Uberhange schlagt sich nun in einer effektiven Energiedichte der Oberflache 
nieder. 

Wir konnen nun den ersten Beitrag zur Hamiltonfunktion aufschreiben. 



r ist die Oberflachensteifigkeit (oder Oberflachenspannung) . Das Integral iiber den Wurzelausdruck 
ist die Flache der Wand. Fiir kleine Gradienten kann man diesen Ausdruck nun nach 
{'Vz{x)y entwickeln. In der niedrigsten Ordnung erhalt man 



Hierbei haben wir den Konstanten Term weggelassen. Dieser beeinfiufit nur die Grundzustandsenergie 
und hat somit kcinen Einflufi auf einen Phaseniibergang. 

Bis hierher bezogen sich unsere Betrachtungen auf ein reines System. Im folgenden wollen 
wir das System derart erweitern, daB wir Unordnung zulassen. Dabei wollen wir zwei 
verschiedene Arten von Unordnung betrachten, Random Bond und Random Field Unordnung. 

Die erste Art der Unordnung kommt dadurch zustande, daB im System an einigen Stellen 
nichtmagnetische Atome sitzen. Im Ising-Modcll modcUicrt man eine solche Unordnung, 
indem man die Kopplungskonstante J durch eine Zufallsvariable ersetzt. 

J ^ Jq + 5J 





(1.1) 
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Abbildung 1.2: Veranschaulichung des Vergrdberungsprozesses: Auf jedem Gitterpunkt sitzt 
ein Spin, der entweder nach oben (+) oder nach unten (— ) zeigt. Die Grenze zwischen 
Bereichen + und — ist die zu untersuchende Wand. In dieser Konfiguration gibt es 
noch Uberhdnge beziehungsweise Einschliisse. Um diese zu beseitigen, mittelt man iiber 
kleinere aber hinreichend grofie Bereiche ( exemplarisch dargestellt durch die gestrichelten 
Ellipsen). Man endet bei einer vergrdberten Beschreibung der Wand, die nun in etwa wie 
die durchgezogene wellige Linie aussieht. 



Jq ist hier eine Konstante, etwa die Kopplungskonstante des reinen Systems. 6 J ist eine 
Zufallsvariable, welche von Ort zu Ort einen anderen Wert besitzt. Ihr Mittelwert ist SJ = 
0. Ubersetzen wir dieses nun in die Formulierung unserer Oberflachenhamiltonfunktion, so 
wird die Oberflachenspannung F eine ZufallsgroBe 

n(x J d^x (r + sr{x, z{x)))y/i + {Vz{x)y. 

Wir entwickeln hier nun wiedcr bis zur ersten Ordnung in ST und (V-z)^. 

n(x J d^x h:{Vz{x)f + 5r{x, z{x)). 

Bei der zweiten Art der Unordnung, der Random Field Unordnung, befindet sich das 
System in einem Zufallsfeld, einem Feld, welches an jedem Ort im Raum einen zufalligen 
Wert annimmt. Dieses bedeutet, daB jeder Spin ein zufalliges 'auBeres' Feld sieht. Im Ising- 
Modell erhalt man hierdurch den zusatzlichen Energiebeitrag HiSi, wobei das Feld 
ist, welches der Spin auf dem i-ten Platz sieht. Die Hamiltonfunktion der Wand erhalt den 
Zusatz 

Hrf - J d'^x J^'^ \z' H{x,z') = J d^x U{x,z{x)). 
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1.1 Ein Oberflachenmodell am Beispiel des Ising-Modells 



Bezeichnen wir nun im Random Bond Fall den Term ST [x, z) mit U {x, z) , so haben wir 
in beiden Fallen dieselbe Gestalt der Hamiltonfunktion. Zusammen mit dem elastischen 
Anteil lautet sie 



Zwischen den beiden Arten der Unordnung unterscheidet man durch die Art des Unordnungskorrelators. 
Die Random Bonds sind voneinander unabhangig, also unkorreliert. Ihr Korrelator ist somit 
lokal. Die Random Fields sind ebenfalls voneinander unabhangig verteilt, woraus auch 
hier ein lokaler Korrelator fiir das Feld H{x, z) folgt. Durch die Integration in z-Richtung 
wird der Korrelator von U{x, z{x)) in dieser Richtung nichtlokal. Er verhalt sich fiir groBe 
Differenzen wie —\z — z'\. 

ZusammengefaBt erhalt man fiir beide Unordnungen 



wobei die Art der Unordnung durch die Funktion R{z) unterschieden wir. Sie hat in den 
beiden Fallen die Gestalt 



Die Unordnung verandert das Verhalten des Systems dramatisch. 1st es im Fall rein Thermischer 
Fluktuationen so, dafi die Energie der Wand nur zu deren Flache proportional ist, so kommt 
nun noch ein Beitrag durch die Unordnung hinzu. Die Energie kann nun abgesenkt werden, 
indem die Wand sich bevorzugt an solchen Stellen befindet, wo die Kopplungskonstante 
zwischen den cinzclncn Spins besonders klcin ist. Der Grundzustand ergibt sich nun aus 
der Konkurrenz der beiden Energieterme, insbesondere mufi er nicht mehr voUstandig eben 
sein, wie im reinen Fall. 

Wir woUen nun in unsere Hamiltonfunktion noch einen dritten Term einfiigen. Durch die 
Kontinuumsbeschreibung des Problems ist eine Eigenschaft des Systems verlorengegangen, 
namlich die Tatsache, dafi die Position der Wand immer zwischen den Spins liegt und somit 
diskretisicrt ist. Wir wollen dem nun Rechnung tragen, indem wir noch ein zusatzliches 
periodisches Potential in ^-Richtung einfiigen. Die Wand mochte bevorzugt in einem Minimum 
dieses Potentials sitzen. Somit fiihren wir hiermit die Diskretisierung quasi durch die 
'Hintertiir' wieder ein. Wir haben also den additiven Beitrag zur Energie 



Vq ist die Amplitude des pcriodischen Potentials, a der 'Gitterabstand' in 2;-Richtung. Hier 
ist es nun angebracht, eine Bemerkung zu den Gitterkonstanten zu machen. Wir haben 
mit dem Ising-Modell auf einem Quadratgitter begonnen. Die Gitterkonstante war dort in 




(1.2) 



U{x,z)U{x',z') =S{x-x')R{z-z'), 



Rrb{z) (X 6{z) 
Rrf{z) oc -\z\ 
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alien Raumrichtungen gleichgroB. Dieses muB nach dem Ubergang zum Kontinuumsmodell 
nicht mehr der Fall sein. Dabei wurden die Fluktuationen der Wand iiber einen gewissen 
Bereich der Wand herausgemittelt. Die Abmessungen dieses Bereiches in den x-Richtungen 
miissen nicht gcnau so groB sein wie die in die z-Riclitung. Wir wollen dahcr zwischcn 
zwei Gitterkonstanten a\\ und a± unterscheiden. Die crstc (a||) soil die Gittcrkonstantc 
parallel zur Oberflache, also die in den 5f-Richtungcn sein. Sie definiert die kurzwellige 
Abschneidefrequenz A — l/a\\ der auftretenden Impulse (Ultraviolett-Cutoff). Die andere 
Gitterkonstante a± ist die Periodenweite des periodischen Potentials. 

Die gesamte Hamiltonfunktion lautet somit 

j d^x ^{yz{x)f - VqCOs + U{x,z{x))^ (1.3) 

Dieses ist die Beschreibung der Oberflache, welche den Betrachtungen dieser Arbeit zugrunde 
hegt. 

1.2 Die Oberflachenrauhigkeit 

Der Hauptgcsichtspunkt dieser Arbeit ist der Rauhigkeitsiibergang von einer glatten zu 
einer rauhen Oberflache. Um dicscn untersuchen zu konnen, ist es zunachst einmal notwendig 
zu definieren, was man unter der Rauhigkeit einer Flache verstehen soil. Eine sinnvoUe 
GroBe hierfiir sind sicherlich die Fluktuationen der Grenzflache um ihren Mittelwert. Die 
Rauhigkeit einer Wand (beziehungsweise genauer gesagt deren Quadrat) wird definiert als 
die mittleren quadratischen Fluktuation der Hohe z{x) (Hohen-Hohen-Korrelation). 

W{\x-x'\) = {[z{x) - z{x')f)^'^ (1.4) 

Wobei (...) fiir die thermische Mittelung und der Querstrich fiir die Unordnungsmittelung 
steht . 

Konvergiert diese Grofie fiir grofie Abstande |x — x'| gegen einen endlichen Wert, so handelt 
es sich um eine glatte Oberflache. Divergiert die GroBe, so ist die Flache rauh. In der rauhen 
Phase verhalt sich die Rauhigkeit fiir groBe L = \x — x'\ wie 

W{L) oc L^. 

Der kritische Exponent C, heiBt Rauhigkeits exponent. 

1.3 Die Behandlung von Systemen mit Unordnung 

Ein zentraler Punkt bei der Behandlung von ungeordneten Systemen ist die korrekte 
Handhabung der Unordnung. Insbesondere muB darauf geachtet werden, daB nur geeignete 
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GroBen gemittelt werden. Diese GroBen heiBen selbstmittelnd. Die Begriindung hierfiir 
liegt darin, daB man zwischen dem wahrscheinlichsten Wert und dem Mittelwert einer 
GroBe unterscheiden muB. Der Mittelwert ist der Wert, den man bei einer Mittelung der 
entsprechenden GroBe erhalt. Fiihrt man eine Messung an einem realen System durch, so 
miBt man jedoch den wahrscheinlichsten Wert. Nun gibt es aber auch ZufallsgroBen, bei 
denen sich der wahrscheinlichste und der mittlere Wert wesentlich unterscheiden. Dieses 
tritt bei stark asymmetrischen Verteilungsfunktionen auf. Die Berechnung des Mittelwertes 
solcher GroBen ist von sogenannten seltenen Ereignissen dominiert. Bei einer Messung 
haben diese aber quasi keine Bedeutung, da sie eben sehr selten sind. 

Eine dieser GroBen ist die Freie Energie 

F = -TlnTr e'^ = -T\nZ. 

Hierbei ist die Boltzmannkonstante, wie auch im folgenden, gleich 1 gesetzt. 

Um diese GroBe iiber die Unordnung zu mitteln, bedient man sich des sogenannten Replikatricks. 
Zunachst stellt man den Logarithmus als einen Grenzwert dar. 



Die n-te Potenz der Zustandssumme Z" wird nun als die Zustandssumme eines n-fach 
vervielfaltigten Systems betrachtet. 



Tr' bedeutet hier, daB die Spur liber die Freiheitsgrade des vervielfaltigten Systems gebildet 
wird. Dieses sind n mal so viele Freiheitsgrade, wie im urspriinglichen System. Tin ist 
die Hamiltonfunktion des gesammten oder replizierten Systems, also aller Teilsysteme. 
Die Teilsysteme, von denen jedes fiir sich dem urspriinglichen System entspricht, sind 
zunachst ungekoppelt. Die Freiheitsgrade der Teilsysteme werden mit der Systemnummer 
des Teilsystems indiziert z^, a = l..n. Die Unordnung hingegen ist unabhangig von der 
Teilsystemnummer oder wie man auch sagt, vom Replikaindex. Diese Zustandsumme des 
replizierten Systems wird nun iiber die (gauBverteilte) Unordnung gemittelt. Hierdurch 
entsteht eine Kopplung zwischen den einzelnen Systemen. 

Problematisch bei der Anwendung dieses Replikatricks ist, daB hierbei die Bildung des 
Grenzwertes lim„^o und der Unordnungsmittelung vertauscht werden. Es ist mathematisch 
nicht sichergestellt, daB man bei beiden Moglichkeiten denselben Grenzwert erhalt. Dieses 
Problem konnen wir hier nicht weiter untersuchen und erwahnen es hier nur der Vollstandigkeit 
halber. 

1.4 Bekannte Resultate 

Das dieser Arbeit zugrundeliegende Problem wurde in dieser Form bereits von Bouchaud 
und Georges mittels eines Variationsverfahrens untersucht Sie approximieren dabei die 



Inx = lim — (x" — 

n^O n 
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ganze replizierte Hamiltonfunktion durch einen quadratischen Ausdruck und untersuchen 
das Problem mit Hilfe der Replikasymmetriebrechung. Fiir Oberflachendimensionen D 
groBer als zwei (Raumdimensionen d groBer als drei) linden sie drei verschiedene Phasen, 
eine glatte, einc glasartig glatte und eine rauhe. In D = 2 Dimensionen linden sie ebenlalls 
drei verschiedene 'Zustande' des Systems. Aul Grund von heuristischen Argumenten, die 
einen Phaseniibergang in 2 Oberflachendimensionen ausschlieBen, argumentieren sie, daB 
die von ihnen gelundenen Grenzlinien Crossoverlinien sein miissen. 

Fiir vereinlachte Modelle gibt es ebenlalls bekannte Resultate. Solche Vereinlachungen 
bestehen darin, daB man entweder das periodische, das ungeordnete oder auch beide 
Potentiale weglaBt. 



1.4.1 Elastische Grenzflache 

Das einlachste System besteht nur aus dem elastischen Anteil der Hamiltonlunktion. 

- j d^'x^ {^z{x)f (1.5) 

Dieses Modell ist exakt losbar, da die Hamiltonlunktion quadratisch im Feld z ist. Die 
Rauhigkeit ergibt sich zu 

Die Impulse q im Integral sind nach unten begrenzt durch 1/L und nach oben durch l/a\\. 
Es gibt hier drei mogliche Falle zu unterscheiden: 



D > 2 In diesem Fall nimmt das Integral liir x, L — > oo einen endlichen Wert an. Die 



Grenzflache ist also nach unserer Deflnition glatt. 



D = 2\ Dieses ist die sogenannte marginale Dimension. Das bedeutet, die groBte Dimension 
bei der die Oberflache liir endliche Temperaturen noch rauh ist. Eine typische Eigenschaft 
der marginalen Dimension ist, daB die Rauhigkeit logarithmisch mit dem Abstand x wachst 
und somit liir x, L — > oo divergiert 

Wl{x) oc In (- 



D < 2 Auch in diesem Fall ist die Oberflache rauh. Im Vergleich zum vorherigen Fall 
wachst die Rauhigkeit allerdings nach einem algebraischen Gesetz mit dem Abstand x. 
Dieses ist also wesenthch schneller als im marginalen Fall. 

W:,{x) oc ^L-^ 
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Im zweiten und dritten Fall erhalt man nur flir T = eine glatte Oberflache. Fiir endliche 
T ist die Wand immer rauh. Fiir den letzten Fall kann man den Rauhigkeitsexponenten ( 
bestimmen. Er betragt 

Cel = (1.6) 



1.4.2 Sinus-Gordon Modell 

Das nachste Modell ist das Sinus-Gordon Problem. Hierbei setzt sich die Hamiltonfunktion 
aus dem elastischen Term und dem periodischen Potential zusammen. 

Hel = jd^'x |^(V^(f))'-K)COs(^^^(f)^| 

Dieses Modell ist nicht mehr exakt zu behandeln. Man muB geeignete Naherungen durchfiihren, 
wie zum Beispiel eine Variationsrechnung oder Renormierungsgruppenrechnung ^ |ll], 



25| . Auch bei diesem Modell ist die marginale Dimension D = 2. Im Fall D > 2 ist die 
Wand somit glatt. Fiir D = 2 andert sich das Verhalten derart, daB es nun einen durch die 
Temperatur getriebenen Phaseniibergang gibt. Unterhalb einer kritischen Temperatur Tc 
ist die Wand glatt, wahrend sie sich oberhalb in einer rauhen Phase befindet. Die Art des 
Phaseniibergangs fallt in die Universalitatsklasse des Kosterlitz Thoules Ubergang. 



1.4.3 Elastische Oberflache in Ungeordneten Systemen 

Ein weiteres bekanntes Modell ist die elastische Grenzflache in einem ungeordneten Medium 
0, beschrieben durch die Hamiltonfunktion 

n = J d^x !^^{Vz{x) f + Uix,z{x)) 

Fiir Dimensionen D < 4 ist dieses System durch einen T = Fixpunkt bestimmt. Das 
System ist durch die Unordnung bedingt selbst bei T = rauh. Der kritische Exponent 
( laBt sich mit Hilfe einer funktionalen Renormierungsgruppenrechnung bestimmen. Man 
erhalt in der 1-loop Approximation: 

A- D 



Crf 3 

Crb ^ (4 - /}) ■ 0.2083 (1.7) 
Der Random Bond Exponent ist ein numerisches Ergebnis. 

In einer Dimension {D = 1) laBt sich das Problem fiir Random Bond Unordnung mit Hilfe 
von Transfermatrizen exakt losen. Man erhalt 

C = | (1.8) 
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was nicht ganz mit dem numerischen Resultat von Fisher iibereinstimmt. Hierbei ist 
zu beachten, daB die Ergebnisse von Fisher im Rahmen der fiihrenden Ordnung einer e 
Entwicklung gewonnen wurden {e = 4 — D). 

1.5 Skalenargumente 

Einen guten Uberbhck erhalt man bereits, wenn man das Skalenverhalten der auftretenden 
GroBen betrachtet. Dabei reskahert man die auftretenden Langen und Impulse derart, 
daB man zu groBeren Langenskalen iibergeht (wir interessieren uns hauptsachhch fiir das 
Verhalten auf groBen Langen). Durch die Forderung der Invarianz der Hamiltonfunktion 
oder die Equilhbierung von Energiebeitragen erhalt man die kritischen Exponenten. 

Als einfiihrendes Beispiel wollen wir das Modell der elastischen Grenzflache ( p. .51 ) betrachten. 
Wir betrachten, wie sich die Energie einer Ausbeulung der Oberflache mit deren Grofie 
andert. Die Langen parallel zur Oberflache skalieren wie L. Somit liefert das Integral 
J d^x einen Beitrag L^. z ist auf dieser Langenskala die Rauhigkeit W{L). Der Gradient 
V liefert eine inverse Langenpotenz L~^. Somit skaliert die gesamte Hamiltonfunktion wie 
T L^~^W^ . Diese Energie muB von der thermischen Energie aufgebracht werden, wir setzen 
den Ausdruck also gleich T und losen nach W auf. 



Wir erhalten denselben Rauhigkeitsexponenten wie in ( |1.6| ) C, = {2 — D)/2. Auch der 
Vorfaktor (T/F)^/^ ist derselbe. 

1.5.1 Das Imry Ma Argument 

Wir wollen nun eine elastische Grenzflache in einem ungeordneten Medium betrachten und 
die Energiebeitrage des elastischen Term und der Unordnung equillibrieren. Dieses fiihrt 
zum Imry Ma Argument. 

Bei diesem Argument betrachtet man die Beitrage zur Energie bei T = 0, die man gewinnt 
und verliert, wenn man einen Bereich mit dem Durchmesser L einer glatten Oberflache um 
eine Distanz W in z-Richtung verformt. Es entsteht dabei ein Energieverlust durch den 
elastischen Term. Dieser ist genau so groB wie im Fall ohne Unordnung, skaliert also wie 



Die Energie, die notwendig ist, um diese Verformung zu realisieren, kann nur durch die 
Unordnung aufgebracht werden. Das ungeordnete Potential verschwindet im Mittel. Wir 
miissen also die Fluktuationen derselben betrachten. Weiterhin muB zwischen Random 
Bond und Random Field Unordnung unterschieden werden. Im ersten Fall gewinnt die 
Oberflache an Energie, indem sie moglichst viele Minima des ungeordneten Potentials 




FL^-W^. 
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erreicht. Die Fluktuationen sind also proportional zur Flache der Ausbeulung und skalieren 
somit wie . 

Der Mechanismus, durch den die Oberfiache im Random Field Fall Energie gewinnen 
kann, funktioniert anders. Gebrochene Bindungen haben nur dann eine Auswirkung auf 
die Energie, wenn sie an der Oberfiache liegen. Im Volumen spiclt es keine RoUe, ob die 
an diescr Bindung liegenden Spins beide nach oben zeigen oder durch die verschobene 
Oberfiache bcidc nach untcn. Die Random Field Unordnung wirkt auf die Einstcllung 
eines jeden Spins und nicht auf die relative Einstellung zweier benachbarter Spins. Wird 
die Oberfiache verriickt, so kommt die Energieanderung dadurch zustande, daB die Spins 
in dcm von der Oberfiache durchwanderten Volumen ihre Einstellung umgedreht haben. 
Die Energiefiuktuationen sind also proportional zu diesem Volumen und skalieren somit 
wie L^W. 

Wir equillibrieren nun die beiden Energiebeitrage. Dabei ist zu beachten, daB der Energiebeitrag 
der Unordnung die Wurzel aus den Fluktuationen ist. Nachdem wir nach der Rauhigkeit W 
aufgelost haben, konnen wir den Rauhigkeitsexponenten fiir beide Arten der Unordnung 
ablesen 

_A-D 

Crb = ^ (1.9) 

Fiir den Random Bond Fall ist die Rauhigkeit geringer als fiir den Random Field Fall. Der 
Grund hierftir ist, daB die Random Bonds nur an der Oberfiache angreifen, die Random 
Fields jedoch am Volumen. 



1.5.2 Das Flory Argument 

Eine weitere Methodc, den Rauhigkeitsexponenten zu bestimmen, ist das Flory Argument. 
Hierbei skaliert man die Langen und fordert, daB die reduzicrtc Hamiltonfunktion Ti/T 
invariant bleibt, sowie daB der elastische und der ungeordnete Anteil gleich skalieren. Durch 
diese Forderungen wird der kritische Exponent C, festgelegt. Die Langen x und z skalieren 
wie X — x/b und z — z/b^. Mit der Tilde ~ werden reskaherte GroBen bezeichnet. Die 
Unordnung skaliert wie die Wurzel des Korrelator. 

C/Rp(x,S) = 6(^-<)/^C/Rp(f,^) 
i/RB(x,5) = 6(^+^)/'C/RB(x,^) 

Somit ergibt sich 

_A-D 
Crf-^- 

A-D 

5 ■ (1.10) 
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Bei Random Field Unordnung stimmen die Exponenten, die man aus Imry Ma und Flory 
Argument gewinnt, iiberein. Anders ist dieses bei Random Bond Unordnung. In einer 
Oberflachendimension haben wir die Moglichkeit, die bei den Rauhigkeitsexponenten Criory = 
3/5 und Cimry Ma = 3/4 mit dem bekannten exakten Resultat Cexakt = 2/3 zu vergleichen. 
Es stellt sich heraus, daB 

Criory ^ Cexakt ^ Cimry Ma- 



Das Imry Ma Argument iiberschatzt den Exponenten, wahrend die Reskalierung der Hamiltonfunktion 

ihn unterschatzt. Eine Begriindung hierfiir zeigen exakte Grundzustandsuntersuchungen 

fiir D = 1 dimensionale Systeme In Abbildung (|1.3|) sind die optimalen Pfade (dieses 

sind D = 1 dimensionale Oberflachen) zwischen einem gemeinsamen Anfangspunkt und 

einem variierenden Endpunkt fiir Random Bond und Random Field Unordnung dargestellt. 

Die Grafiken zeigen einen markanten Unterschied zwischen Random Field und Random 

Bond Unordnung. Im Random Field Fall "teilen" sich benachbarte Pfade wesentlich friiher 

als im Random Bond Fall. Benachbarte Pfade enthalten bei Random Bond Unordnung 

viele gemeinsame Segmente. Daher muB man die Annahme iiber den Energiegewinn durch 

die Unordnung, welche dem Imry Ma Argument zugrunde liegt, als eine obere Grenze 

betrachten. Der Rauhigkeitsexponent Cimry Ma ist damit ebenfalls eine obere Schranke. 
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(a) 
(b) 

Abbildung 1.3: Die optimalen Grundzustdnde von Benachbarten Grenzfldchen bei Random 
Field (a) und Random Bond (b) Unordnung in einer Oberfldchendimension: Es werden 
Fldchen betrachtet, die einen gemeinsamen Anfangspunkt besitzen (Spitze des Dreiecks) und 
deren Endpunkte auf der Basis des Dreiecks liegen. Deutlich zusehen ist, dafi im Random 
Field Fall benachbarte Fldchen cine wesentlich geringere Anzahl gemeinsamer Segmente 



haben als im Random Bond Fall. (Abbildung aus 



In diesem Kapitel wird die Oberflache mit Hilfe der funktionalen Renormierung untersucht. 
Sie unterscheidet sich von der "normalen" beziehungsweise, genauer gesagt, von der Parameter 
Renormierung dadurch, dafi nicht nur die Kopplungskonstanten in der Hamiltonfunktion 
renormiert werden, sondern dafi sich die ganze funktionale Form der Beitrage zur Hamiltonfunktion 
andert. Auf der formalen Seite haben wir bei der Funktionalen Renormierung anstatt von 
gewohnlichen Differentialgleichungen fiir die Kopplungskonstanten, wie bei der Parameter 
Renormierung, nun partielle Differentialgleichungen fiir die funktionalen Anteile der Hamiltonfunktion, 
wie beispielsweise des periodischen Potentials. 

Die funktionale Renormierungsgruppentransformation besteht aus drei Teilschritten. Zuerst 
wird das Feld z{x) in einen langsam und einen schnellfluktuierenden Anteil aufgespaltet. 
Dieses entspricht der Zerlegung des Impulsraumes in eine Kugelschale mit den Impulsen 
A/6 < Ig*! < A und in eine Kugel mit \q\ < A/6 zerlegt. A ist hierbei der obere Cutoff fiir die 
Impulse, der durch das zugrundeliegende Gitter zustande kommt: A oc l/a||. 6 > 1 ist ein 
dimensionsloser Parameter. Uber die Impulse in der Kugelschale wird integriert. Die dabei 
entstehenden Terme mufi man zu geeigneten Grofien in der verbleibenden Hamiltonfunktion, 
die nur noch von den 'kleinen' Impulsen abhangt, addieren. Dadurch entstehen effektive 
Grofien. Man erhalt so aus der urspriinglichen Hamiltonfunktion mit einem Cutoff von A 
eine, bei der sich der Cutoff auf A/6 reduziert hat. Als dritten Schritt skaliert man alle 
Langen und Impulse so, dafi der urspriingliche Cutoff wiederhergestellt wird. 

Wir betrachten eine D-dimensionale Grenzflache in einem D+l-dimensionalen Raum. Die 
Flache wird durch einen D-dimensionalen Vektor x parametrisiert. Die Auslenkung der 
Flache sei durch das einkomponentige Feld z{x) gegeben. Wir beschreiben die Flache mit 
der folgenden reduzierten Hamiltonfunktion: 

n = j d''x^\h:{Vz{x) f -Vocos{^z{x)^ (1.11) 

Reduzierte Hamiltonfunktion bedeutet, dafi diese durch die Temperatur T dividiert wurde. 
Im weiteren Verlauf dieses Kapitels soil der Begriff Hamiltonfunktion immer die reduzierte 
Hamiltonfunktion meinen. 
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Der erste Term beschreibt die elastische Energie der Oberflache, der zweite ist ein periodisches 
Potential, welches die diskrete Struktur des Gitters simuliert. Der letzte Term ist das 
Unordnungspotential, welches gauBverteilt ist und durch die beiden Momente 



V{z,x) = 



V{z,x)V{z',x') = S{x - x')R{z - z') 

charakterisiert wird. Alle hoheren Momente lasscn sich durch das Wick Theorem auf diese 
beiden zuriickfiihren. Der Querstrich iiber den GroBen symbohsiert hierbei die Unordnungsmittelung. 
Mit der Funktion R{z) wird die Art der Unordnung modelliert. Fiir die beiden in dieser 
Arbeit betrachteten Unordnungformen hat sie die Gestalt 



R{z) oc 



— \z\ random field 
6{z) random bond 



Die Unordnungsmittelung behandeln wir mit Hilfe des Replikatricks. Die replizierte Hamiltonfunktion 
hat folgende Form 

= Jd^xf^Y.^ ^{Vz"{x)f-Vocos(^^z-{X)^ 

-^E^(^"(^)-A^))| (1-12) 

Die griechischen Buchstaben a und P sind die Replikaindizes, sie laufen von 1 bis n. Im 
Weiteren wird die Abkiirzung V{z) :— V^) cos (27r2;/aj_) verwendet. 



1.6 Skalierung 

Zunachst woUen wir uns einen Uberblick verschaffen und das Verhalten der einzelnen 
Anteile in der Hamiltonfunktion bei Reskaherung der auftretenden Langen betrachten. 
Wir ersetzen die GroBen x und z durch x und z mit 

X — e'x, z — e^^z. 

Die Auslenkung z der Flache muB mit dem Rauhigkeitsexponenten C, skaliert werden, da 
sie wie die Rauhigkeit skahert. Hier und im folgenden deutet die Tilde ~ an, daB es sich um 
reskalierte GroBen handelt. Diese setzen wir in die replizierte Hamiltonfunktion ein. Wir 
fordern nun, dafi die Kopplungskonstante F des elastischen Terms unverandert bleibt und 
erhalten somit das Skalenverhalten der Temperatur. 

/ d^xe™-^ ( ^ze^^-A^ = f d^x^ ( ^zX 
J 2r \dx J J 2T \dx J 
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Mit diesem T bestimmen wir nun ganz analog R{z) und V{z) 

Die hieraus folgenden Flufigleichungen lauten 

D-2C)f (1.12.a) 
D -AC)R + CzR' (1.12.b) 

2C)V^ + aV^' (1.12.C) 

An diesen FluBgleichungen kann man bereits zwei wesentliche Dinge ablesen. Die Temperatur 
T flieBt fiir D > 2 — 2^ gegen null. In dem fiir uns interessanten Bereich D >2 haben wir 
es also mit einem T = Fixpunkt zu tun. Wir miissen also bei dcr folgenden Funktionalen 
Renormierung nur solche Terme mitnehmen, die auch im Grenzwert T einen Beitrag 
liefern, die anderen fallen weg. An der FluBgleichung (1.12.b) fiir R sieht man, daB die obere 
kritische Dimension 4 ist. D.h. oberhalb von 4 Dimensionen ist die Unordnung irrelevant. 
Die Grenzflache bleibt also wie fiir D > 2 Dimensionen im reinen Fall glatt. 



dT 



dl ^ 

dl ^ 



1.7 Funktionale Renormierungsgruppe 

Wir beginnen die Funktionale Renormierung in der iiblichen Art und Weise mit der Zerlegung 
des Feldes z{x) in den lang- und kurzwelligen Anteil. 



2(f) = 

z'^{x) = 


z<(x) 

1 


q 


VVd 


z'^{x) = 


1 


Y.^ ^q 

q 


Wd 



ist dabei einc Abkiirzung fiir die Summe iiber Impulse mit A/b < \q\ < A, Ylq f^r die 
iiber Impulse mit \q\ < A/b. 

Der elastische Term ist diagonal in diesen beiden Teilfeldern. Das periodische Potential 
sowie der Unordungsanteil miissen noch nach diesen Feldern zerlegt werden. Beide Terme 
werden nach den kurzwelligen Feldern bis zur zweiten Ordnung entwickelt. 

V{z) ^V{z< + z>) = V {z<) + V (^<) ^> + ]-V" {z<) + 0(^>3) 
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Der Strich ' symbolisiert die Ableitungen nach z. Bei der Funktion R ist noch darauf zu 
achten, dafi sie nicht diagonal in den Replikaindizes ist. Zur Vereinfachung fiihren wir die 
Abkiirzung 2;"^ :— z"^ — ein. 



Bei der Summation iiber die Replikaindizes kann man die Symmetrie der Funktion K 

R{z) = R{-z) ^ R'{z) = -R'{-z) => R"{z) = R"{-z) 
sowie z"'^ — —z^°' ausnutzen. Man erhalt somit 

\ 7 / 



Wir fiihren folgend Fouriertransformationen ein: 







V"'^{q) := 






Jr. 




[ d^'x 



Vol 
18 



^^e^^^^^V"{z-<{x)) 

f d^'x e'^ R' (z"^<(a;)) 

J VVol 



1.7 Funktionale Renormierungsgruppe 



und erhalten die Zerlegung der Hamiltonfunktion 
mit 

a y q q,q' y 

-2^E ( E'^'"'(^1-I- E'^"'(^"'+^"'>i4j 

a,/3 y q q,q' J 



Ti,{z°'^'^) hat dieselbe funktionale Gestalt wie die urspriingliche Hamiltonfunktion. Der 
einzige Unterschied zwischen bcidcn ist, daB diese nur noch von den kleinen Impulsen 
abhangt. Hq[z"^) ist eine quadratische Hamiltonfunktion, die nur groBe Impulse enthalt 
und zudem diagonal in den Rcplikaindizes ist. Tii hangt von alien erlaubten 

Impulsen ab. Somit werden hierdurch die Anteile mit kleinen und groBen Impulsen gekoppelt. 
Durch die Renormierungsgruppentransformation entstehen aus diesem Teil der Hamiltonfunktion 
die Korrekturen zu Tilz"^^). 

Hiermit konnen wir nun die kurzwelligen Freiheitsgrade ausintegrieren, das heiBt wir bilden 
die Zustandssumme, wobei wir nur iiber die Feldanteile mit den groBen Impulsen integrieren. 
Das Resultat interpretieren wir als exp(— 7ii.enormiert(^"^^))- 



Y J-oo \ /Hoiz>) 

= e-^(^")Zo>cxp(e-«i(^"'^")-l 



no{z>),c 



Um die Kumulanten zu berechnen, benotigt man die Erwartungswerte der verschiedenen 
Potenzen von z^. Da die Hamiltonfunktion Tio quadratisch und diagonal in den Impulsen 
q ist, kann man das Wick Theorem anwenden. Die Erwartungswerte mit einer ungeraden 
Anzahl von z verschwinden. Die mit einer geraden werden auf den quadratischen Erwartungwert 
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zuruckgefiihrt, wobei 



ist. 



a 13 



Ho(z>) 



T 



In erster Ordnung erhalt man: 



T 



22-2 



q,q' a,f3 



rg2 



= E'E(iv-(0) + ^Q-(0))^ 



q a. 



Die hier auftretende Summe iiber die Impulse q nahern wir durch ein Integral: 



> 1 



Vol 



d^'q 1 



Vol K 



.A 

D / dq (f' 

JKIh 



= Vol Kn 



'A/6 
1 



2-D\ 



D-2 



Kd ist dabei die Oberflache dcr D-dimensionalen Einheitskugel dividiert durch (27r)^. 
Dieser Term zeigt fiir D > 2 und 6 — > oo keine Divergenz. Somit konnen die Beitrage der 
ersten Ordnung vernachlassigt werden. 

Als nachstes wird die zweite Ordnung betrachtet: 



a y q q,qi J 

E I E' fi'"'(<?)^l - E' '?')44 



a,/3 



9,9 



Ho{z>),C 



Durch Ausmultiplizicrcn des Quadrats [. . . ]^ entstehen 10 verschicdcnc Summandcn, die 
als Kandidaten fiir Korrckturen in Frage kommen. Vier sind von der Ordnung und fallen 
somit bei der Mittelung weg. Drei sind quadratisch in z: 



\-2 E E^ ^'"(^"^^'"'(^"")^"'"'^?.-«-' 



22-2 



T 

IV 



2T 



OL q 
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°':I3 q,q' ^ 



a/3 q-q-' 

a' 



a,/3 q ^ 

Bei diesen drei Beitragen taucht wieder 1/g^ auf, so dafi sie keine Divergenz fiir 6 — > oo 
zeigen. Sie werden also vernachlassigt. 

Der Term, der quadratisch in V" ist, ist von der Ordnung T°. Wie wir bereits weiter oben 
gesehen haben, handelt es sich bei diesem Problem um einen T = Fixpunkt. Die Terme, 
die einen Beitrag liefern, miissen also von derselben Ordnung in T sein, wie die GroBen 
in der urspriinglichen replizierten Hamiltonfunktion. Der Term des periodischen Potentials 
ist von der Ordnung . Somit wird der in V" quadratische Term vernachlassigt. 

Der quadratische Term in Q sieht folgendermaBen aus: 



/3' 



= ir^E E'«°'('?+'?')«°'(-'?-'?')^i 

Als erstes ersetzen wir die Summe iiber q und q' durch entsprechende Integrale und fiihren 
k — q + q' ein: 

rri q^ q^ 



r {k + qf 

Der Kringel ° am Integral zeigt an, daB k so gewahlt sein soil, dafi q und q' in der 
Kugelschale mit A/6 < |g | < A liegen. Wir nehmen an, daB die Q's nur einen Beitrag 
fiir kleine k liefern. Dann nahern wir l/[k + q)'^ ~ ^/(f- Somit vereinfacht sich das Integral 
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zu 



Vol^ 




ci^xQ^^Q"^— A-^(6^ - 1). 



Dieser Term zeigt eine Divergenz fiir D < A und h ^ oo. Wir setzen nun die Q"^ ein: 



Hier haben wir drei verschiedene Beitrage. AUe drei sind von der Ordnung T~^. Die ersten 
beiden enthalten zwei Replikaindizes, somit liefern sie einen Beitrag zur Renormierung von 
R. Der letzte Term hangt von drei Replikaindizes ab. Er wiirde also eine 'Drei-Replika- 
Kopplung', die in der urspriinglichen Hamiltonfunktion nicht vorhanden ist, renormieren. 
Es bleibt die Frage, ob dieser Term eine drei Replika Kopplung generieren kann. Wir 
betrachten dazu, wie eine Kopplung zwischen drei Replikas in der replizierten Hamiltonfunktion 
zustande kommen konnte. Die Mischterme entstehen durch die Unordnungsmittelung. Um 
einen Term mit drei Replikaindizes zu generieren, miiBte die Unordnungsverteilung nocli 
nicht gauBsclie Anteile enthalten. AuBer dem quadratischen miiBte es noch andere Momente 
der Unordnungsverteilung geben, die sich nicht durch das quadratisch Moment ausdriicken 
lassen. Eine hieraus entstehende Kopplung von drei Replikas ware aber von der Ordnung 
T~^0. Somit ist dieser Term fiir einen T = Fixpunkt irrelevant, da er von der Ordnung 
ist. Folglich vernachlassigen wir ihn. 

Es bleibt noch ein Mischterm zwischen dem ungeordneten und periodischen Potential zu 
betrachten. 




7 



1 



T T 



2T3 



Tq^ Tq 



,'2 




Hierbei wurde wieder k = q + q' eingefiihrt und die Annahme gemacht, daB das Integral 
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1.7 Funktionale Renormierungsgruppe 



nur Beitrage fiir kleine k liefert. In der Ortsdarstellung haben wir nun 



e 



Der erste Summand hangt von einem Replikaindex ab. Er ist von der Ordnung T~^. 
Somit handelt es sich um eine Korrektur zu V. Der zweite enthalt eine Summe iiber zwei 
Replikaindizes, er konnte also einen Beitrag zu R liefern. Dazu hat er allerdings die falsche 
Ordnung in T. Er fallt also weg. 

Damit haben wir alle Terme in dieser Ordnung betrachtet. Wir erhalten insgesamt die 
effektiven GroBen 



R"'(z) 



2R"{z)R"{0) 



Veff = V{z) - -^-A-^{b^ -l)V"{z)R"{0). 



:i.l5) 
:i.l6) 



Fiir die Temperatur T ergeben sich in dieser Ordnung der Entwicklung keine Korrekturen. 



Die replizierte Hamiltonfunktion hat nun wieder die Gestalt von ( |1.12| ), wobei das periodische 
Potential und der Unordnungskorrelator durch die effektiven GroBen ersetzt worden sind. 
Der obere Cutoff hat sich von A auf A/6 reduziert. 



Wir reskalieren die auftretenden Langen wie in Abschnitt |1.6| und stellen somit den urspriinglichen 
Cutoff wieder her. Mit einem infinitesimal gewahlten / = ln(6) erhalten wir folgendes 
FluBgleichungssystem: 



dT 

'm 

dR 
'dl 
dV 
'dl 



— = (2 - D - 2C)T 



— = {t-K)R+CzB! 



KnA 
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R \z) - 2R"{z)R"{0) 



(2 - 20V + (zV - -^A-^V"{~z)R"{0) 



(i.ie.a) 
;i.i6.b) 

(1.16.C) 



Die ersten beiden Gleichungen entsprechen genau denen von D.S. Fisher 0. Neu hinzugekommen 
ist die dritte Gleichung fiir das periodische Potential. Eine erweiterte Form dieser Gleichung 



taucht bei J. Kierfeld beim Problem gekoppelter Schichtsupraleiter auf [|13|. Wir definieren 
nun noch 

D{1) := !^m, 
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1 Das Modell 



um uns von den physikalisch unbedeutenden Konstanten in den FluBgleichungen zu befreien 
und erhalten 



dT 

— = {2-D- 2C)T 
^ = (2 - 2C)V + CzV' - ^V"{~z)D"{0). 



D \~z)-2D"{~z)D"{Q) 



(1.17.a) 
;i.l7.b) 
(1.17.C) 



1.8 Fixpunkte der Flufigleichungen 

Um die verschiedenen Phasen zu bestimmen, ist es notig, die Fixpunkte der FluBgleichungen 
zu bestimmen. Besonders einfach ist dieses ftir die Temperatur. T = ist hier ein Fixpunkt. 
Da die zugehorige FluBgleichung eine lineare, gewohnliche Differentialgleichung erster Ordnung 
ist, kann man sie auch voUstandig losen 

f{l) = T{1 = 0)e(2-^-2O-'. 

Somit wird die Temperatur fiir groBe I und Dimensionen D > 2—2(, was groBen Langenskalen 
entspricht, gegen null getrieben. T = entspricht einem stabilen Fixpunkt. Er beschreibt 
somit eine Phase und keinen Phaseniibergang. Dieses Resultat haben wir schon bei der 
Herleitung der FluBgleichung fiir das periodische Potential und den Unordnungskorrelator 
verwendet . 

Schwieriger ist es bei den anderen beiden FluBgleichungen, da es sich hier um partielle 
Differentialgleichungen handelt. Weiterhin miissen wir noch anstelle eines Anfangswertes, 
wie bei der Temperatur, Randbedingungen, die durch die urspriinglichen, unrenormierten 
GroBen gegeben sind, beriicksichtigen. Fiir das periodische Potential V ist die Randbedingung, 
daB die Funktion periodisch ist. Die Randbedingung des Unordnungskorrelator ist bei den 
beiden Unordnungsarten unterschiedlich. Im Random Field Fall ist die Randbedingung die 
Asymptotik des Korrelators fiir groBe z, R{z) oc —\z\. Im Random Bond Fall hingegen die 
Tatsache, daB die Korrelationen kurzreichweitig sind. 

Wir betrachten zunachst die FluBgleichung fiir den Unordnungskorrelator. Dieses ist genau 
dieselbe FluBgleichung, die D.S. Fisher in seiner Arbeit zu einem ahnlichen Problem, 
namlich dem ohne Periodischen Potential, gefunden hat . Wir erhalten also fiir die beiden 
Unordnungsarten dieselben Fixpunkte wie er. 

Zunachst betrachten wir den Random Field Fall. Setzen wir hier das asymptotische Verhalten 
des Unordnungskorrelators ein, so ist damit der Rauhigkeitsexponent ( bereits festgelegt. 

C = ^ (1.19) 
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1.8 Fixpunkte der Flufigleichungen 



Mit diesem ( haben wir nun die Gleichung 

~Di~z) + '-~zD'{~z) + \ [b"\z) - 2b"{z)b"{Q)\ = 

zu losen. Leiten wir diese Gleichung nun einmal nach z ab und dividieren sie danach durch 
D"{z), so erhalten wir eine elementar integrable Gleichung. Diese ausgefiihrt ergibt nun 
eine implizite Gleichung fiir die zweite Ableitung der Fixpunktfunktion D* (z) 

A* In (^-^-^^ +^* + D*"iz) = (1.20) 

wobci A* := —D*"{0) gcsctzt wurde. A* ist durch die obige Gleichung noch nicht bestimmt 
und parametrisiert somit eine ganze Familie von Fixpunktfunktionen, die sich durch einen 
unterschiedlichen LangenmaBstab in 2;-Richtung unterscheiden. Eine natiirliche Wahl fiir 
A* ist 

A* = e. (1.21) 

Im Falle von Random Bond Unordnung kann man die Fixpunktfunktion nur numerisch 
bestimmen. Fiir grofie z verhalt sich hier der Korrelator am Fixpunkt wie 

b*{z) o^ez~^-^h-^. (1.22) 
Die numerische Berechnung liefert den Rauhigkeitsexponenten 

C ^ 0.2083e (1.23) 

Als letztes miissen wir nun die Fixpunktfunktion fiir das periodische Potential finden. Die 
Unordnung geht hier durch den Faktor D"{0) ein. Setzen wir hier V = 0,so stellen wir fest, 
daB dieses ein Fixpunkt der FluBgleichung ist. Wir miissen jetzt noch die Stabilitatseigenschaften 
dieses Fixpunktes untersuchen. Die Reskalierung in z-Richtung bringt den schwierigen 
Term (zV'{z) in die FluBgleichung hinein. Um diesen zu vermeiden, betrachten wir das 
zwar renormierte aber nicht reskalierte effektive Potential. 

Veff{z) = V{z)-l-V"{z)D"{0){¥-l) 
Mit Qo :— A/b erhalten wir somit die FluBgleichung 

Da es sich hierbei um GroBen handelt, die nicht reskaliert wurden, der gefundene Fixpunkt 
D aber eine reskalierte GroBe ist, muB die Skalierung noch riickgangig gemacht werden. 

b{z) = 6^-^^D(6^5) 
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1 Das Modell 



Weiterhin sehen wir, daB sich der Cosinus des periodischen Potentials durch die zweite 
Ableitung von V nach z reproduziert. Wir nehmen nun an, dai3 durch die Renormierung 
nur der Vorfaktor Vq des periodischen Potentials beeinfluBt wird. Somit erhalten wir hier 
eine FluBgleichung fiir eben diesen Vorfaktor Vq: 



Wir erhalten also eine lineare Differentialgleichung der Ordnung 1. Diese konnen wir nun 
durch Trennung der Variablen losen. Wie integrieren die resultierende Gleichung qo = 1/L 
bis qo = A und erhalten 



Wir sehen, daB das periodische Potential auf groBen Langenskalen immer irrelevant wird. 
Das bedeutet, daB unsere Grenzflache immer rauh ist. Die hier vorgestellte Rechnung 
licfcrt also kcincn Phascniibcrgang zwischen einer glatten und cincr rauhen Phase. Die 
Art der vorhandenen Unordnung geht hier nur iiber den Rauhigkeitsexponenten ( ein. 
Somit verschwindet das periodische Potential fiir den Random Field Fall schneller als fiir 
den Random Bond Fall. Dieses vertragt sich mit der Vorstellung, daB die Random Field 
Unordnung eine starkere Unordnung als die Random Bond Unordnung ist. Bei der Random 
Bond Unordnung spielen nur die Defekte in einer relativ kleinen Umgebung der Grenzflache 
eine Rolle. Die zusatzlichen Energiebeitrage der Zufallsfelder erhalt man, indem man die 
Flache von z = ±00 in das System hineinschiebt. Es handelt sich hierbei also um eine Art 
Volumeneffekt. Dies spiegelt sich auch in der Korrelationsfunktion wieder, fiir Random 
Fields ist sie langreichweitig, fiir Random Bonds hingegen lokal. 




Vo{L) = Vo{a\\) ■ cxp 




(1.24) 
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Kapitel 2 

Variationsrechnung 



In diesem Kapitel wollen wir den EinfluB des periodischen Potentials mit Hilfe der Variationsrechnung 
untersuclien. Bei diesem Verfaliren wird die urspriingliclie Hamiltonfunktion durcli eine 
einfachere ersetzt. Diese vereinfaclite Hamiltonfunktion enthalt nocli Parameter, die zu 
bestimmen sind. Dieses geschieht durcli die Minimierung eines geeigneten Variations- 
Funktionals (Anhang 

:= Fo + {H- no)n, ■ 

Hierbei ist Ti. die urspriingliclie und Tio die Variations-Hamiltonfunktion. Fq ist die Freie 
Energie der Variations-Hamiltonfunktion: 

Fa = -TlnTr e"^ 



Auch liier reclinen wir wieder in e = 4 — D Dimensionen. Wir ersetzen das periodische 
Potential durch ein quadratisches, im folgenden Masseterm genannt. Der zu variierende 
Parameter ist die Kopplungskonstante dieses Terms, die Masse m. Wir liaben also 



n 

Ho 



d^'x <J ^ iVz{x)f - Vo cos (^z{x)^ + U{z{x), x) 

r 1 

d^x \ - i\/z{x)f + -m^z^{x) + U{z{x),x) 



(2.1) 
(2.2) 



Hierbei handelt es sich wieder um replizierte Hamiltonfunktionen. Den Index n haben wir 
hier und im folgenden weggelassen. 

Wir wollen nun das Variations-Funktional minimieren. Dazu berechnen wir 



d 



^ T^i-^lii^o + ^Tr 
om om I Zo 



{n - no)e 



T 
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2 Variationsrechnung 



wobei = Tr exp(— Tio/T") die Zustandssumme bezliglich der Hamiltonfunktion Tio ist. 
Diesen Ausdruck gleich null. Wir benutzen die Tatsache, daB nur Hq und Zq von m 
abhangen: 







Zn Zf 



+ — Tr 



•'0 

■9^0 

dm 



dZ. 



dm 



l + -{n-Ho)]e 



'Hp 
T 



Wir berechnen weiter 

dZp _ d 

dm dm 

und erhalten somit 



■Hp 

Tr e T 



Tr 



1 dHo 
T dm 



■e T 



Zq /dHo 

T \ dm 







dm 



■Ho 



dHo 

dm 



Ho 



Wo 



Ho 



T 



9m 



Wo 



Ho 



Ho 



Diese konnen wir noch als Kumulanten-Erwartungswert zusammenfassen und erhalten 
somit nach der Multiplikation mit T die Selbstkonsistenzgleichung 







{Ti - Ho 



dm 



(2.3) 



Ho,C 



Wir miissen nun diese Kumulante berechnen. Wir nehmen dazu an, daB das Feld z{x) 
durch die Variations-Hamiltonfunktion gauBverteilt ist. Dieses gilt nach D.S. Fisher ^ 
zumindestens fiir T = 0. Durch diese Annahme konnen wir das Wick-Theorem verwenden 
(Anhang [T^). 

/27r 







d^x — Vo cos ( — z(x] 



-m^z^ix) 



- I d^x d^y 



TC 



VqCos I — z{x) H — m z (x) 



d^y mz'^{ 
mz'^{y) 



Ho,C 



Ho,C 



d X d y m {z{x)z{y))^^ {m -Vq\ 



/27ry^-^((.(.))^),^ 



Hinreichend fiir die Erfiillung dieser Gleichung ist das Verschwinden des Ausdrucks in 
der geschweiften Klammer. Hierbei ist zu beachten, daB eine Funktion von m ist. 
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2 Variationsrechnung 



Wir wollen hier natiirlich wieder ein unordnungsgemitteltes Ergebnis haben. Daher miissen 
wir hier den unordnungsgemittelten Erwartungswert von verwenden. Somit lautet die 
Selbstkonsistenzgleichung 

= Vo i—j e "-0. (2.4) 



Wir berechnen {z'^{x))^^ mit Hilfe der replizierten Hamiltonfunktion. 



Hierzu entwickeln wir die replizierte Hamiltonfunktion bis zur zweiten Ordnung in z". 
Dieses liefert unter der Annahme, daB das z gauBverteilt ist, ein korrektes Ergebnis. Des 
weiteren transformieren wir noch die Hamiltonfunktion in die Fourierdarstellung, da diese 
diagonal in den Impulsen g ist. 

^''-f = / (0 S:^ ((r-?^ + '"^)^'' - ^-|^fi"(«)") ;1^^r 

=:(G-i(9))"'' 

Die Mittelung besteht nun im wesentlichen in der Invertierung der Matrix G^^ (Anhang 
|1.3| ). Wir erhalten 

Der erste Summand ist dabei der thermische Beitrag zur Rauliigkeit. Aus dem reinen 
Modell wissen wir, daB der oberlialb von zwei Oberflachendimensionen eine verschwindende 
Rauhigkeit liefert. Der zweite Summand ist der Unordnungsanteil an der Rauliigkeit. Sein 
Beitrag ist von der Ordnung T^. Er liefert also auch bei T = einen Beitrag. 

So wie bis jetzt formuliert, wird die Unordnung "nur" storungstheoretisch behandelt. Dieses 
fiihrt, wie wir bereits aus dem Fall ohne periodisches Potential wissen, zu ungenauen 
Resultaten aufgrund von vielen lokalen Minima der Hamiltonfunktion. Wir wollen daher 
fiir die Unordnung die Ergebnisse der Renormierungsgruppe benutzen. Dabei wird die 
Kopplungskonstante durch eine effektive Kopplungskonstante ersetzt. Diese effektive Kopplungskonstante 
ist die GroBe, die in der Renormierungsgruppe nach Ausintegration der schnellen Freiheitsgrade 
definiert wurde. Sie ist dadurch skalenabhangig geworden, weshalb wir besser von einem 
Kopplungsterm sprechen. 
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2 Variationsrechnung 



2.1 Drei- Parameter- Approximation 



Eine besonders schone Methode, den effektiven Kopplungsterm zu berechnen, besteht in der 

Hierbei wird 



Drei-Parameter- Approximation von T. Nattermann und H. Leschhorn [16, 2J 
die funktionale Renormierung auf eine bestimmte Familie von Funktionen eingeschrankt. 
Die so gefundenen FluBgleichungen flir die Parameter dieser Funktionenfamilie sind elementar 
integrierbar, so daB wir nicht nur den Fixpunkt der Renormierungsgruppe erhalten, sondern 
auch den FluB der Parameter dorthin. Ein Nachteil dieser Methode soil nicht verschwiegen 
werden. Sie ist nur auf solche Unordnungskorrelatoren anwendbar, deren Integral iiber die 
ganze z-Achse nicht divergiert. Dieses ist fiir das Random Bond Problem gewahrleistet, 
wenn man den anfanglichen Unordnungskorrelator durch eine GauBfunktion darstellt. 

Wir beginnen mit einer geeigneten Definition von D{z) 



Diz) 



A 



r{u) 





\u\ 




\u\ 



(2.6) 



Mit dieser Definition fiihren wir zwei Parameter A und ^ ein, die renormiert werden sollen. ^ 
ist die Korrelationslange der Unordnung. A gibt die Starke der Unordnung an. Die Funktion 
r ist kurzreichweitig. Wir verwenden hier spater fiir Rechnungen eine GauBfunktion. Diesen 
Ansatz fiir D setzen wir nun in die FluBgleichung (l.lT.a.b) fiir D aus der funktionalen 
Renormierung ein. Wir erhalten 



A fz 



Az 



J ^ 



di ^ 

z\ A2 



r — 



2 r" ( - ) r"(0) 



Von dieser Gleichung bestimmt man nun zwei verschiedene Momente beziiglich z. Die 
einzelnen Summanden der Gleichung sind gerade Funktionen von z, da der Unordnungskorrelator 
eine gerade Funktion ist. Die Momente, die wir auswahlen, miissen gerade sein, da die 
ungeraden verschwinden. Aus diesen beiden Gleichungen fiir die Momente konnen wir nun 
die FluBgleichungen fiir die Parameter A und ( gewinnen. Wir erhalten 



[e - 5C + CAg]A 



dA 
'dl 



(2.7) 



Hier tauchen drei neue GroBen auf: ca, und g. Die ersten beiden bestimmt man aus der 
folgenden Gleichung 



Ca + mc^ 



J du u"'r"{u)[r"{u) - 2r"(0)] 



{21 
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2.1 Drei-Parameter- Approximation 



Dabei setzt man fiir m ganze gerade Zahl ein, die entsprechenden Momente. g ist eine 
dimensionslose Kopplungskonstante aus der Kombination der beiden Parameter A und ^: 
g = A/^^. Aus den Gleichungen ( f^.T] ) laBt sich eine FluBgleichung fiir g herleiten: 

^ = eg + [cA- 5cd (2.9) 
Der Fixpunkt dieser Gleichung berechnet sich zu 

9* = ^ • (2.10) 

Man kann die FluBgleichung nun von / = bis / = In(AL) integrieren 

9(0 I ri 
dg y r- = / dl' 



und erhalt 



9. eg(l-j,) JO 



~g(l) = M!^!! ^ (2.11) 

g*+go{e^'-l) ^ ' 



mit go = g{l = 0) = ^ 
so 



Damit konnen wir nun auch die FluBgleichungen fiir A und ^ losen: 



Ao 



A io I g* + go{e''-l) 



Co 



^ Jo I 'g* + go{e''-l] 



Wir erhalten 



A(AL) = AoG'^(L) 

aAL)=^oG^^{L), (2.12) 
wobei wir zur Vereinfachung c := c^/c^ und 

G{L) := 1 + ^ [(A^)^ - 1] 

definiert haben. 

Wir ersetzen nun die Kopplungskonstante —R"{0) in ( |2.5| ) durch die aus dem Ansatz ( |2.6| ) 
herriihrende. Die GroBen A und ^ werden durch die unskalierten Versionen der effektiven 
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2 Variationsrechnung 



GroBen ( p.l2| ) ersetzt. Wir konnen nun {z"^) fiir T = berechnen: 

'l\ 1 



Ao 
so 



l/L 



(Fg^ -|- 777,2)2 ^ 



(2.13) 



i^'i) ist hier wiederum das Volumen der D-dimensionalen Einheitskugel im Impulsraum, 
dividiert durch (27r)^. Die Masse m bewirkt nun, dafi die kleinen Impulse abgeschnitten 
werden. Wir haben also in einem unendlich groBen System einen unteren Cutoff von m/ \/V. 

"A 




- 1 



Somit konnen wir nun die Selbstkonsistenzgleichung aufstellen. Es ist im weiteren niitzlich, 
die konkurrierenden physikalischen Beitrage durch miteinander vergleiclibare Langenskalen 
auszudriicken. 

Zuerst einmal haben wir bereits drei Langen in der Gleichung stehen, a±, und 1/A . a± 
ist die Periode des gittersimulierenden Potentials. Es ist somit eine Lange in z-Richtung. 
A ist der obere Impulscutoff, somit ist 1/A = a\\ die Gitterkonstante in x-Richtung. 

.^0 ist das Maximum aus der Breite der Grenzflache (bei einer Blochwand ist das zum 
Beispiel die Langenskala, auf der sich die Spins um 180° drehen) und der Korrelationslange 
der Unordnung. Auch dieses ist eine Lange in 2;-Richtung. Diese Lange gibt an, wie stark 
man die Wand verschieben kann, ohne daB sich der Energiebeitrag der Unordnung drastisch 
andert. Nehmen wir zum Beispiel an, wir haben eine Wand der Breite I, die an einem 
Punktdefekt, welcher zunachst in der Mitte liegen soil, eine gewisse Energie AE gewinnt. 
Wir konnen nun diese Wand um ±d/2 verschieben. Der Punktdefekt bleibt dadurch in der 
Wand, wodurch sich der Energiegewinn nicht wesentlich andert. Analog ist die Betrachtung 
einer diinnen Wand, welche mit Defektstellen endlicher Breite wechselwirkt. (Abbildung 



Die nachste auftretende Lange in unserem Problem beschreibt die Unordnung. Es ist 



die sogenannte Larkin-Lange [|1^]. Das ist die Langenskala, ab der die Unordnung in 
der Art relevant wird, daB man sie nicht mehr mit der Storungstheorie behandeln kann. 
Dieses ist dann der Fall, wenn die Rauhigkeit W der Wand von der GroBenordnung der 
Gitterkonstante an wird. 



W^{L) 



2Ao 

He 



^1 -1 



Setzen wir nun Qo und g* (x e ein, so erhalten wir 



^29o I L 



^0 _ 
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2.1 Drei-Parameter- Approximation 




Abbildung 2.1: Die Korrelationsldnge C,o- Es gibt zwei Extremfdlle. Nehmen wir zundchst (a) 
an, die Punktdefekte hahen eine gewisse Breite (ii^ der Abbildung durch die Gaufikurve 
dargestellt) und die Position der Wand sei sehr scharf definiert. In diesem Fall kann die 
Position der Wand im Intervall C,o variieren ohne das sich der Energieanteil der Unordnung 
dramatisch dndert. Wir haben in der Skizze drei mdgliche Positionen als gestrichelte Linien 
eingezeichnet. Im anderen Fall (b) sind Punktdefekte sehr schmal und die Wand hat eine 
Breite ^o- Auch hier ist es mdglich, die Wand in diesem Intervall zu verschieben. In der 
Abbildung (b) ist die durchgezogene Linie die Position der Defektstelle. Die gestrichelten 
Linien zeigen exemplarisch zwei mdgliche Positionen der Wand. 



Wir definieren daher die Larkin-Lange als 

n* , \ \ 

(2.14) 




Wir sehen, daB eine schwache Unordnung einer groBen Larkin-Lange und eine starke 
einer kleinen Larkin-Lange entspricht. Weiterhin wird die Larkinlange mit wachsender 
Korrelationslange kleiner. Die Auswirkungen der Unordnung sind also um so starker, je 
breiter die Defektstellen sind. 

Die letzte Langenskala ist dadurch gegeben, daB die Wand eine gewisse Distanz in x- 
Richtung benotigt, um von einem Minimum des periodischen Potentials zum nachsten 
zu gelangen [Q. Diese Lange ist gepragt von der Konkurrenz zwischen der elastischen 
Energie und dem periodischen Potential. Der elastische Term versucht die Wand zu glatten, 
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2 Variationsrechnung 



Z 



Minimum des periodisclien Potentiais 



/ 



Maximum des periodisclien Potentials 



Wand 





a 




W 



X 



Abbildung 2.2: Die Ldngenverhdltnisse, die beim Wechsel der Oberfldche von einem 
Minimum zum ndchsten des Periodischen Potentials auftreten. Eine solche Stelle wird als 
Stufe oder Kink bezeichnet. 

wodurch der Ubergang zwischen zwei Minima des periodischen Potentials in die Lange 
gezogen wird. Das periodische Potential wirkt diesem entgegen. Der Energieverlust durch 
dieses Potential ist proportional zur Breite dieser Stufe. Die typische Breite der Stufe ist 
nun dadurch bestimmt, dafi die Energiebeitrage der beiden Terme von der gleichen GroBe 
sind. Der Beitrag der elastischen Energie ist proportional zur Breite multipliziert mit 
dem typischen Gradienten zum Quadrat (a^/w)^. Mit den Vorfaktoren Vq und F erhalten 
wir also 



Zuletzt skalieren wir noch den Variationsparameter m zu einer dimensionslosen GroBe p 
um. Wir dividieren dazu m durch den Vorfaktor der Exponentialfunktion auf der rechten 



Der physikalisch relevante Bereich, in dem p liegen kann, liegt zwischen und 1. Bei 
p = 1 entspricht die Variations-Hamiltonfunktion Tig einer Entwicklung des periodischen 
Potentials in der ursprlinglichen Hamiltonfunktion Ti bis zur quadratischen Ordnung in 
z. Werte von p, die groBer als 1 sind, wiirden bedeuten, daB das periodische Potential 
die Rauhigkeit unserer Wand starker unterdriickt als das Parabelpotential der Variations- 
Hamiltonfunktion. Dieses ist physikalisch unsinnig: Der Energieverlust pro Wandflache fiir 




(2.15) 



Seite von Gleichung ( |2.4| ) und definieren dieses als p^ . 

a_|_ 1 
p := m ■= 




(2.16) 
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2.2 Auswertung der Selbstkonsistenzgleichung 



eine rauhe Wand bleibt im periodischen Potential immer endlich. Im Parabelpotential wird 
er hingegen flir beliebig groBe Auslenkungen beliebig groB. 

Wir haben nun alle GroBen, um die Selbstkonsistenzgleichung (|2.4|) aufzuschreiben: 




(2.17) 

Die auftretenden GroBen sind nun entweder dimensionslos, wie p und die Konstanten c und 
cj, oder Langenverhaltnisse. Die Langen werden in den natiirlichen Einheiten des Systems 
gemessen. Dieses sind die Gitterkonstanten in der entsprechenden Richtungen: Ll und w 
werden mit a\\, mit a± gemessen. 

Wir miissen jetzt noch die Konstanten und c mit Hilfe von Gleichung ( p.8|) bestimmen. 
Um diese exakt zu bestimmen, ware es notwendig, die exakte Fixpunktfunktion fiir r{u) 
einzusetzen. Diese ist jedoch analytisch nicht bestimmbar und liegt somit nur numerisch 
vor. Wir miissen also eine geeignete Naherung fiir diese Funktion finden. Hieraus folgt, daB 
die Konstanten nun von der Wahl der Momente m abhangt. Wir setzen an: 




Der Nachteil dieser Funktion ist, daB sie im Gegensatz zur tatsachlichen Fixpunktfunktion 
am Ursprung analytisch ist. Die richtige Fixpunktfunktion zeigt hier eine Divergenz in 
der vierten Ableitung. Man sollte daher m 7^ wahlen. Die Abhangigkeit von m bei der 
genaherten Fixpunktfunktion ist allerdings sehr gering. Mit den Momenten 2 und 4 erhalten 
wir 

c = 0.1885... 
= 2.1105... 

Der hieraus resultierende Rauhigkeitsexponent 

C = ^ 0.2078e 

5 — c 

kommt recht nahe an Fishers Ergebnis [0 von 0.2083e. 



2.2 Auswertung der Selbstkonsistenzgleichung 

Wir woUen nun die Selbstkonsistenzgleichung ( |2.17| ) analysieren und deren Losungen finden. 
Zunachst stellen wir ein paar allgemeine Uberlegungen an. Wir betrachten die urspriingliche 
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Form der Selbstkonsistenzgleichung ( p.4| ). Mit dem reskalierten Variationsparameter p 
lautet sie 

= exp I ^{z^{x)) 

Weil ihr Argument immer negativ ist, ist die Exponentialfunktion ihrerseits immer kleiner 
als eins. Somit ist die Bedingung, daB fiir eine physikalische sinnvolle Losung p kleiner 
als eins sein muB, immer erfiillt. Bei der Berechnung von {z{x)) wurde eine Naherung 
gemaclit, die niclit allgemein giiltig ist. Im Integral iiber die Impulse (|2.13| ) ist die untere 
Integralgrenze durch die Masse modifiziert worden. Dieses ist jedoch nur dann sinnvoU, 
wenn m/\/T kleiner als die obere Abschneidefrequenz der Impulse ist. Fiir groBere Werte 
wird die Naherung extrem schlecht und bewirkt zum Beispiel, daB das Argument der 
Exponentialfunktion in der Selbstkonsistenzgleichung ( |2.17| ) negativ werden kann, wodurch 



Losungen p > 1 folgen. Die Bedingung m/^/T < A iiber setzt sich in die Bedingung 

fiir die neu definierten GroBen. Bei Gleichheit ist der Wert der Parameter erreicht, bei dem 
die Exponentialfunktion in Gleichung ( |2.17| ) gleich eins wird. 



Zunachst woUen wir den Fall D > 4 betrachten. e ist in diesem Fall negativ, so daB 
die Exponentialfunktion fiir alle p analytisch ist. Sie fallt mit p aufgrund des negativen 
Vorzeichens im Exponenten und konvergiert gegen null fiir p oo. Die Gleichung hat 
somit also immer eine Losung p > 0. Ein Beispiel fiir konkrete Werte der vorkommenden 
GroBen zeigt Abbildung ^.3| . Die Grenzflache ist fiir diesen Fall also immer glatt, wie wir 
erwartet haben. 

Kommen wir nun zum Fall D < 4. In diesem Fall ist e groBer als null. Das Argument der 
Exponentialfunktion ist an der Stelle p = nicht analytisch. Es divergiert gegen —oo fiir 
p —>■ 0, wodurch die Exponentialfunktion selber null wird. p = ist also immer eine Losung 
des Systems. Es verschwinden alle Ableitungen der Funktion im Grenzwert p — 0. Es kann 
also nicht passieren, daB eine weitere Losung aus der Null "herauswandert" , wie es zum 
Beispiel bei der Mean-Field Behandlung des reinen Ising-Systems der Fall ist. 

Fiir groBerwerdende p hat man das Problem, das irgendwann die Bedingung (|2.18 ) nicht 



mehr erfiillt ist. In diesem Fall ist die Naherung zu verwerfen und man muB das Integral 
numerisch behandeln. 



Das typische Aussehen der Selbstkonsistenzgleichung zeigt Abbildung |2]J. Dort ist ein 
Fall gezeigt, wo die Gleichung drei Losungen hat, eine Losung fiir p = und zwei fiir 
nicht verschwindende p. Die p = Losung existiert fiir alle moglichen Werte der Parameter 
Ll, ^ und w. Die Existenz der anderen beiden Werte pi < p2 ist von diesen Parametern 
abhangig. Die Wahl der Parameter bestimmt wie "schnell" bzw. steil die Exponentialfunktion 
wachst. 

Wir woUen nun untersuchen, wie sich die Losungen pi und p2 bei kleinen Variation der 
Parameter Ll, und w verandern. Erhohen wir w, dieses entspricht einem schwacheren 
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Abbildung 2.3: Die Ahhildung zeigt das typische Aussehen der Selbstkonsistenzgleichung 
fiir Dimensionen D > 4. Hier wurde speziell D = 5, w/a\\ = 4, Li/a^ = 2 und 
^q/0'± = 5 gewdhlt. Auf der Abszisse wird , auf der Ordinate p^ beziehungsweise die 
Exponentialfunktion aufgetragen. Der Schnittpunkt der beiden Kurven ist die Losung der 
Gleichung. 

periodischen Potential, so wird p2 kleiner und pi groBer. Wir schlieBen daraus, daB die 
groBere der beiden Losungen, also p2, die pliysikalische Losung ist: Ein kleineres p unterdriickt 
die Rauhigkeit weniger als ein groBes. Weiterhin bewirkt ein schwacheres periodisclies 
Potential eine grofiere Rauhigkeit als ein starkes. Dieses beiden Tatsachen treffen fiir p2 zu. 
Analog ist das Verhalten bei Anderung von und L^. Erlioht man oder senkt man L^, 
so wird in beiden Fallen p2 kleiner und pi groBer. Dieses paBt ebenfalls in das pliysikalische 
Bild: Eine starkere Unordnung soUte eine groBere Rauhigkeit und somit ein kleineres p 
zur Folge haben. Bei beiden Fallen erhoht man die Unordnung, einmal, indem man ihren 
"EinfluBbereich" vergroBert (^o erhoht) und zum anderen indem man den numerischen 
Wert der Unordnung erhoht {Ll kleiner macht). 

Der Phaseniibergang ist genau die Wertekombination von ^, Li und w bei der eine Losung 
p 7^ zum ersten mal auftritt. Betrachten wie die Gleichung ( p.4|) , so stellen wir fest, daB, 
wenn die Exponentialfunktion fiir irgendeinen Wert der Masse m ungleich null ist, wir mit 
einem geniigend groBen Vq einen Schnittpunkt zwischen dieser und p^ erreichen konnen. Die 
Frage nach der Existenz eines Phaseniibergangs reduziert sich hierdurch auf die Frage, unter 
welchen Bedingungen die Exponentialfunktion fiir alle p den Wert null annimmt. Dieses 
ist wiederum genau dann der Fall, wenn das Argument der Exponentialfunktion nach —oo 
divergiert. Die Divergenz muB wegen dem oben gesagten unabhangig fiir alle p auftreten. 
Dieses ist der Fall, wenn die Larkin-Lange gegen null geht oder die Korrelationslange 
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Abbildung 2.4: Hier ist eine Situation fiir D < 4 gezeigt. Die Auftragungsart ist dieselbe wie 
in Abbildung Hier hat sich das Verhalten der Exponentialfunktion drastisch gedndert, 
so dafi es bis zu drei Schnittpunkten gibt. 



^0 divergiert. Beides entspricht einer unendlich groBen Unordnung und ist somit nicht 
realisierbar. Eine andere Moglichkeit ware eine Divergenz der Lange w. Man kann dieses 
durch Vo — > oder F — > cxd erreichen. Es handelt sich also auch hier um nicht allzu sinnvoUe 
Grenzwerte, da zum einen kein periodisches Potential vorhanden ist und zum anderen die 
Wand unendlich steif und somit immer glatt ist. 

Beschranken wir uns auf endliche, von null verschiedene Werte der Parameter, so ist die 
Exponentialfunktion auBer bei p = von null verschieden. Somit gibt es fiir diese Wahl der 
Parameter immer einen durch das periodische Potential V getriebenen Phaseniibergang. 

Wir stellen nun die Lange w, die die Starke des periodischen Potentials beschreibt, als 
Funktion der anderen beiden Langen und Ll dar. Fiir fest vorgegebene Werte der 
Korrelations- und Larkin-Lange bestimmen wir den maximalen Wert w{^o,Ll) bei dem 



es in Abbildung noch einen Schnittpunkt zwischen den eingezeichneten Kurven gibt. In 
diesem Fall tangiert die Exponentialfunktion die Gerade. Dieses maximale w stellt somit 
einen Phaseniibergang zwischen einer glatten Phase, bei der p ^ ist, und einer rauhen 
mit p = dar. Der Phaseniibergang ist diskontinuierlich, das heiBt, p springt von einem 
endlichen Wert auf null. 



Hieraus ergibt sich ein Phasendiagramm, wie es Abbildung |2.5| zeigt. Die dargestellte Flache 
ist die Grenzflache zwischen der rauhen (oberhalb) und der glatten Phase (unterhalb). 
Wir sehen, daB eine starkere Unordnung, das heiBt eine kleinere Larkin-Lange Ll, eine 
starkeres periodisches Potential benotigt (kleineres w), um in eine glatte Phase zu gelangen. 
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Abbildung 2.5: Das Phasendiagramm: unterhalh der dargestellten Fldche ist der 
Ordnungparameter p 0. Dieses bedeutet, dafi sich das System in der glatten Phase 
befindet. Oberhalb der eingezeichneten Fldche ist das System rauh. Bei der Berechnung 
dieses Diagramms wurde das Integral i\2.1S ) numerisch ausgewertet. 



Bemerkenswert ist die Abhangigkeit des Phaseniibergangs von der Korrelationslange ^o- 
Ist diese grofi gegeniiber der Gitterkonstanten des periodischen Potentials a_L, so ist die 
Abhangigkeit sehr gering. Fiir Werte ^ < a± zeigt sich demgegeniiber eine sehr starke 
Abhangigkeit. In diesem Fall reicht schon ein verhaltnismaBig kleines Potential, um die 
Wand zu glatten. Ist die Korrelationslangen von der GroBenordnung des Gitterkonstante 
oder groBer, so scheint die Unordnung einen Teil des periodischen Potentials "herauszumitteln" , 
weshalb der numerische Wert Vq groBer sein muB, um eine glatte Wand zu erhalten. 

Bei Random Field Unordnung haben wir die Moglichkeit, das Fixpunktergebniss der Unordnung 



zur Auswertung von Gleichung p.4| zu verwenden. Wir erhalten eine FluBgleichung, die 



keinen Parameter enthalt, welcher die Starke der Unordnung genauer charakterisiert. 



= Vq{ — 1 ■ exp 

Die Uberlegungen bei Random Bond Unordnung beziiglich des periodischen Potentials sind 
weiterhin giiltig. Somit haben wir auch hier einen durch Vq getriebenen Phaseniibergang. 



27r\^ 
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Fiir ein Vq, welches kleiner als ein kritischen Vc ist, ist die Oberflache rauh, fiir ein Vq groBer 
als Vc glatt. Dieser Phaseniibergang ist ebenfalls diskontinuierlich. 

Wir woUen nun versuchen, ob man auch aus der Renormierungsgruppenrechnung einen 
Phaseniibergang herleiten kann. Wir andern dazu die rephzierte Hamihonfunktion ( |1.12| ) 
derart ab, daB wir noch einen Masseterm J d^x l/2m'^z^{x) einfiigen. 



K a 



^ (V^"(f))' + ]-m\z'\x)f - ^ocos i—z'^ix] 
2 z \0'i. 



(2.19) 



Mit dieser neuen Hamihonfunktion fiihren wir nun die funktionale Renormierung genau wie 
in Kapitel |1.5.2| durch. Das Einzige, was sich andert, ist der zu verwendende Propagator. 



Wir ersetzen in Gleichung ( |1.14| ) die beiden Propagatoren und fahren mit denselben Naherungen 
fort. Letztendhch bleibt ein Integral der Art 



2^2 



iibrig. Dieses Integral ist uns schon aus der Variationsrechnung bekannt. Es bewirkt, daB die 
untere Impulsabschneidefrequenz sich bei hinreichend groBer Masse m auf m/\/T erhoht. 
Die FluBgleichungen 16.4 b) und ( |1.16.a| .c) andern sich hierdurch der Art, daB die Terme, 
die die zweiten Ableitungen von V beziehungsweise R enthalten, in diesem Fall wegfallen. 



f = (e-4C)i^ + a^' + l^A- 



2R!'{~z)R!'{Q) 



(2.20) 



G ist die Stufenfunktion, die fiir positive Argumente den Wert 1 annimmt und sonst null 
ist. Sie bewirkt hier die Fallunterscheidung zwischen massivem und masselosen Fall. Fiir 
/ — s> 00 fallt der letzte Term bei einer endlichen Masse immer weg. 



Was bedeutet das nun fiir das System? Ohne die Masse m liefert das Integral ( |1.14[ ) im 

thermodynamischen Limes fiir den Grenzwert I — > 00 die divergenten Beitrage, welche 

die Unordnung als auch das peridische Potential auf nichttriviale Art renormieren. Diese 

Divergenzen treten mit Masse nicht mehr auf, da die Impulse nun nach unten beschrankt 

sind. Das bedeutet, daB die Renormierung bei einem endlichen Wert des FluBparameters I 

stoppt. Dieses / entspricht einer endlichen Langenskala, ab der sich die effektiven Kopplungskonstanten 

des Systems nicht mehr andern. Die Oberflache bleibt in diesem Fall also glatt. 
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Die Frage, die noch zu klaren bleibt, ist, welchen Wert hat die Masse m? Der Term, der 
uns diesen leifern kann, ist das periodische Potential. In einer einfachen Storungsrechnung 
wiirde uns die zweite Ableitung an der Stelle z = diesen Term liefern 



g2 



V{z) 



/2n\' 



= Vo{- 

z=0 \(^±J. 



Wir verwenden hier, um die Masse zu bestimmen, ein mit dem massiven Propagator 



renormiertes Vq. Das Vorgehen ist dasselbe wie in Abschnitt 1.8 mit dem Unterschied 



daB die Masse mitgenommen wird. Somit enden wir bei der FluBgleichung flir Vq 



In der Stufenfunktion wurde = A/go ersetzt. Die folgende Impulsintegration wird wieder 
durch die Masse nach unten begrenzt. Somit erhalten wir 



27rV / /27rV 



e 



™:^ = m^j -expl-^^j i^l 1^1 -111 (222) 




Mit diese Gleichung ist nun die Masse m selbstkonsistent zu bestimmen. Es ist die Geleichung, 
welche wir erhalten wiirden, wenn wir im vorangehenden Kapitel den unordnungsgemittelten 
Erwartungswert (z^) nicht mit der drei Parameter Approximation, sondern 'nur' mit dem 
Fixpunktresultat der Renormierungsgruppe ftir die Unordnung, berechnet hatten. Was 
diese Methode nicht liefert, ist der Crossover von kleinen zu groBen Langenskalen. Diesen 
erhalt man nur dann, wenn man die FluBgleichungen der Funktionalen Renormierung 
simultan integriert. Das ist aber nicht ohne weitere Vereinfachungen, wie eben der drei 
Parameter Approximation, moglich. Somit hat sich der Kreis geschlossen. 
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Kapitel 3 
Zusammenfassung 



In dieser Arbeit wurde eine elastische Oberflache in einem ungeordneten Medium mit 
einem zusatzlichen periodischen Potential untersucht. Ziel dieser Arbeit war es, die Existenz 
eines Phaseniibergangs zwischen einer rauhen und einer glatten Phase zu beleuchten. Die 
obere kritische Dimension fiir dieses System ist D=4, genau wie im Fall ohne periodisches 
Potential . Oberhalb dieser Dimension ist sowohl das periodische Potential als auch die 
Unordnung irrelevant. 

Zunachst wurde das Problem mithilfe der funktionalen Renormierung untersucht. Dabei 
wird die Abweichung von der kritischen Dimension als kleiner Parameter e = 4 ~ D 
genommen, und nach diesem bis zur fiihrenden Ordnung entwickelt. Hieraus resultieren 
FluBgleichungen fiir die Temperatur T, das periodische Potential und die Unordnung. Das 
System wird durch einen T = Fixpunkt beschrieben. Die FluBgleichungen fiir die anderen 
beiden GroBen sind partielle Differentialgleichungen. Die Beschreibung der Unordnung wird 
nicht durch das periodische Potential beeinfluBt. Somit ergibt sich hierfiir derselbe Fixpunkt 
wie bei D.S. Fisher |^. Die FluBgleichung des periodischen Potentials laBt sich auf einer 
FluBgleichung der Kopplungskonstante Vq dieses Potentials reduzieren. Diese wiederum 
liefert, eine stabilen Fixpunkt Vq = 0. Somit ist das periodische Potential immer irrelevant. 

Eine zweite Methode, mit der das Problem behandelt wurde, ist die Variationsmethode. Das 
periodische Potential wurde durch ein quadratisches ersetzt. Die dazugehorige Kopplungskonstante 
m wurde dabei so bestimmt, daB ein geeignetes Variations- Funktional minimal wird. Die 
Unordnung wurde bei dieser Rechnung weiterhin mithilfe der Renormierungsgruppe behandelt. 
Diese wurde derart abgeandert, daB der Funktionenraum, in dem sich der Korrelator der 
Unordnung befindet, derart beschrankt wurde, daB sich die Renormierungsgruppe auf die 
Renormierung dreier Konstanten reduziert. Die drei Konstanten sind die Temperatur T, 
die starke der Unordnung A und die Breite der Defektstellen ^. Diese Approximation bietet 
sowohl Vor- als auch Nachteile: Der wichtigste Vorteil ist, daB nun der Ubergang (Crossover) 
von kleinen zu groBen Langenskalen berechnet werden kann. Der groBte Nachteil ist, daB 
hiermit nur lokale Unordnungskorrelatoren, wie sie bei der Random Bond Unordnung 
auftreten, betrachtet werden konnen. Das Ergebnis dieser Rechnungen ist, daB es Unterhalb 
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von D = A Oberflachendimensionen einen Phaseniibergang erster Art zwischen einer 
rauhen und einer glatten Phase gibt. Was diese Rechnung nicht liefert, ist eine untere 
kritische Dimension, ab der das periodische Potential immer irrelevant ist. 

In einem weiteren Abschnitt wurde versucht, auch aus der Renormierungsgruppenrechnung 
einen Phaseniibergang herzuleiten. Zu diesem Zweck ist die Hamiltonfunktion durch einen 
Masseterm erweitert worden. Dieser Term bewirkt, dafi die Renormierung auf einer durch 
die Masse bestimmten Langenskala zum erliegen kommt. Versucht man nun diese Masse aus 
den im urspriinglichen Modell vorhanden GroBen zu extrahieren, so endet man wiederum 
bei einer selbstkonsistent zu losenden Gleichung. Im Gegensatz zur Variationsrechnung 
enthalt diese nur die Fixpunktwerte der Unordnung und keinen Crossover. 
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Anhang 



1.1 Das Variations-Funktional 

Angenommen, wir haben zwei Hamiltonfunktionen Ti und Hq- Wir benotigen ein Funktional 
dieser beiden Funktionen, welches fiir Ti = Tio sein Minimum annimmt. 

Durch die beiden Hamiltonfunktionen werden zwei Verteilungsfunktionen definiert: 



P = — 7^ — Z = Tr e 



-H/T 



Po 

Wir betrachten nun 



g— Ho/T 



Po • (Inpo - Inp) = Po • In — 

P . 

Weiterhin benutzen wir die bekannte Ungleichung fiir Logarithmen 

\nx < X — 1 =^ lnx>l 

X 

und erhalten 

Po • (Inpo - Inp) > Po - p. (1.1) 

Die Spur iiber die rechte Seite dieser Ungleichung verschwindet, da beides auf eins normierte 
Verteilungsfunktionen sind. Somit ist die Spur iiber die linke Seite immer grofier gleich null. 



< Tr{po ■ (Inpo - Inp)} 
= Iv{po-(-lnZo-^ + lnZ+^)} 

^lnZ-]nZo + ^{n-no)no 
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(. . . bedeutet, daB der Erwartungswert mit Tio gebildet wird. Multiplizieren wir dieses 
Funktional mit T und nehmen nur die Anteile mit, die von Tio abhangen, so erhalten wir 
ein flir die Variationsrechnung geeignetes Funktional 



wobei Fq = —T\nZo die Freie Energie des durch Ho beschriebenen Systems ist. Die 
Ungleichung (|1.1| ) ist genau dann eine Gleichung, wenn p = po ist. Somit ist das Variations- 
Funktional genau dann minimal, wenn Ti = TYq gilt- Ebenfalls sind in diesem Fall die drei 
Freien Energien F™, Fq und F gleich. 



1.2 Erwartungswert von cos{27vz/a 



Wir woUen den Erwartungswert von cos(27r2;/ a) beziiglich einer in z quadratischen Hamiltonfunktion 
7i berechnen. 

/27r \\ Tr cos(27rz/a)e-^/^ 
cos — z ' ^ — 



Der dazugehorige Boltzmannfaktor hat eine gauBsche Form. Somit ist das Wick Theorem 
anwendbar. Zu Beginn zerlegen wir den Kosinus in seine Reihenentwicklung 



COS I — z 

a 



n 



Der Erwartungswert von {z'^"')^ berechnet sich nun mit dem Wick Theorem: Man muB die 
Erwartungswerte aller moglichen Paarungen von z berechnen. Da alle z gleich sind ist dieses 
besonders einfach, namlich gleich (-z^)" multipliziert mit der Anzahl der Moglichkeiten, die 
z zu paaren. Fiir das erste z hat man 2n — 1 mogliche Partner. Es bleiben 2n — 2 z iibrig. 
Somit gibt es fiir das zweite Paar 2n — 3 Moglichkeiten. Fiihrt man dieses Verfahren fiir 
alle z durch, so erhalt man insgesammt {2n — 1)!! Moglichkeiten. 



cos I — z 

a 



H 



2\" 
H 



Zuletzt wollen wir versuchen, durch geschickte Umformungen die Summe wieder aufzusummieren. 
Wir zerlegen die Fakultat in zwei Doppelfakultaten. Die ungerade der beiden kiirzen wir 
gegen die bereits vorhandene Doppelfakultat. Die gerade zerlegen wir in 2"n!. 
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cos I — z 

a 



n 
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1.3 Matrixinversion 



Diese Summe identifizieren wir als Exponentialfunktion. 




Analog berechnet sich der Erwartungswert 



2 /^27r 

Zi COS Z2 

\ a 



rtfi 



Weil hier eine Kumulate berechnet wird, ist zu beachten, daB zunachst die beiben zi mit 
jeweils einem Z2 gepaart werden. Man erhalt: 



27r 

Zl cos I Z2 

a 



n,c 



1 f2n 

a 



n,c 



1.3 Matrixinversion 

Gegeben sei eine quadratische n x n Matrix 

M"^ = a5"^ + b, 



[1.2) 



von der die inverse Matrix zu bestimmen ist. Wir nehmen an, daB die Inverse dieselbe 
Struktur besitzt. 



:i.3) 



Es muB gelten 

In = M ■ = M-^ ■ M, 
wobei In die n x n Einheitsmatrix ist. Fiir die einzelnen Elemente der Matrix gilt daher 



Nach Einsetzten von (|1.2|) und (|1.3| ) erhalten wir 

= ^[a ■ X ■ d""^ ■ 5^^ + a ■ y ■ 5'''< + h ■ X ■ 5^^ + h ■ y] 

7 

= a ■ X ■ 5"^ + a- y + b- x + n- b- y. 
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Dieses Gleichungssystem hat die eindeutige Losung 

1 

X — — 
a 

b 

^ a • {a + n ■ b) 



Somit lautet die gefundene Inverse 



a a ■ [a + n ■ 
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